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Abstract 

An improved version of a previously validated linearized Euler equation solver is used to 
compute the noise generated by coannular supersonic jets. Results for a single supersonic 
jet are compared to the results from both a normal velocity profile and an inverted velocity 
profile supersonic jet. 

Introduction 

The full, compressible Navier-Stokes equations govern the process of sound generation 
and propagation to the far field. To solve these equations without resorting to modeling the 
turbulent quantities. Direct Numerical Simulation (DNS) may be employed. However, the 
resolution required to accurately simulate high-Reynolds-number turbulent flows make 
DNS calculations impractical due to current computer limitations. 

Therefore, Mankbadi et. al. 1 - 2 proposed the extension of the Large Eddy Simulation (LES) 
approach for use in the prediction of sound generation and propagation. In this approach. 
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the Navier-Stokes equations are filtered into large-scale components, which are calculated 
directly, and small-scale components, which are modeled. In a LES calculation, the noise 
generated by the small-scale components is not captured; however, it is known that the 
large-scale structures are more efficient than the small ones at radiating noise (e.g., Refs. 3- 
7). The LES approach has been used successfully to predict the noise radiated by a 
supersonic jet, but it is still a CPU-intensive approach, especially for three-dimensional 
calculations of far-field noise. However, this limitation may be bypassed by using either a 
Linearized Euler Equation (LEE) solver 8 or a Kirchhoff method 9 ' 12 to calculate the far-field 
noise radiation. 

The present work is concerned with exploring the use of the less computationally 
demanding LEE solver for supersonic jet noise predictions. The LEE approach neglects 
both viscosity and nonlinear effects. The viscous effects can be neglected since the large- 
scale dynamics in free shear flows are essentially inviscid (e.g.. Ref. 13). Nonlinearity, 
however, seems to be important (e.g., Ref. 14). 

Yet, much of the physics can be obtained by considering the linear equations. Several 
attempts have succeeded in studying the physics of jet noise based on a simplified form of 
the Euler equations (e.g., Ref. 15-18). The linearized Euler equations simultaneously 
describes the near field, where the sound is generated, and the propagation of sound to the 
far field. As such, the problem of matching the near field hydrodynamic disturbances to 
the far field acoustic propagation does not arise. The linearized Euler equations fully 
account for non-parallel flow effects and for the simultaneous presence of non-discrete 
frequencies. 

In the past, the linearized Euler equations have been used to compute the instability noise 
generated by supersonic jets 19 - 20 . Parametric studies have been performed to investigate 
the effects of jet temperature 21 , and the effect of random inflow disturbances has been 
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shown 22 . Other researchers have added nonlinear terms to the linearized Euler equations, 
but their magnitude and effect are uncertain 23 . 

Governing Equations 

Starting from the full Navier-Stokes equations in conservative form, neglecting viscosity, 
linearizing about a given mean flow 
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and separating the azimuthal modes, the linearized Euler equations may be written in 
cylindrical coordinates as: 


where n is the azimuthal mode and: 
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Velocities are normalized by the reference jet exit centerline velocity U c , time by D/U c , 
density by the jet exit centerline value p e , and pressure by p e U c “. Here, D is the reference 
jet nozzle diameter. 
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Reference 

Normal Velocity 
Profile 

Inverted Velocity 
Profile 

U inner (ni/s) 

1330.0 

1534.6 

823.7 

Uouter (m/s) 


1227.7 

1441.4 

Tinner (°K) 

1100.0 

1692.3 

1678.4 

ToutcrCK) 


846.2 

808.1 

Mj nn er 

2.0 

1.9 

1.0 

M 0 uter 


2.1 

2.5 


Table 1 : Supersonic coannular jet test cases 


In Eq. (2), a periodic distribution is assumed for the azimuthal direction; for example, 
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In this formulation, each additional mode calculated adds another set of equations to be 
solved. The advantages to this formulation are: lowered storage and computation, 
improved centerline behavior, and improved boundary condition specification. 
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Reference 
(inner -> ») 

Normal Velocity 
Profile 
(inner -> °°) 

Inverted Velocity 
Profile 
(outer -> 00 ) 

Mck-H 

1.97 

(6 = 59.6°) 

2.03 

(0 = 60.5°) 

2.19 

(0 = 62.8°) 

Mc^supersonic 

1.32 

(0 = 40.6°) 

1.33 

(0 = 41.3°) 

1.56 

(0 = 50.3°) 

Mcl subsonic 

0.659 

0.631 

0.939 


Table 2: Predicted Instabilities and Approximate Emission Angles 

Mean Flow 

The mean flows used in these test cases were computed using a Parabolic Navier-Stokes 
(PNS) numerical solver provided by Dr. Milo Dahl, which is documented in Ref. 24. 
There were three test cases computed, shown in Table 1. 

In these test cases, the thrust and mass flow of the three jets are held constant. The 
coannular jets have an area ratio A 2 /A 1 of 1.25, and a radius ratio x^x\ of 1.5. The Normal 
Velocity Profile test case corresponds to Dahl’s case 5, and the Inverted Velocity Profile 
test case corresponds to Dahl's case 24. 

In heated supersonic jets, there are three possible types of instability waves: the classical 
Kelvin-Helmholtz instability waves, supersonic instability waves, and subsonic instability 
waves. The convective Mach numbers for the three families of instabilities are given in 
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Seiner, et. al.'s 25 work as: 
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When the convective Mach number is greater than unity, that instability may be sustained 
by the jet. An estimate of the Mach wave emission angle is given by: 

•“"•’fe) (12) 

Table 2 gives a listing of the convective Mach numbers for these three families of waves as 
well as the estimated Mach wave emission angle. Notice that the shear layer between the 
inner and outer cores of the coannular jets may also support instabilities. 

Velocity profiles for these three cases are given in Figure 1, and the velocity magnitude 
along the radial center of each jet are given in Figure 2. 

Numerical Algorithm 

The code uses a new MacCormack-type solver, which is formally fourth order accurate in 
time and space. This extension of the MacCormack scheme uses Tam and Webb's 
Dispersion Relation Preserving methodology in its derivation. This scheme is described in 
greater detail in Ref. 26, and has been validated against the 2-4 scheme, which is another 
MacCormack-type solver. The 2-4 scheme was described by Gottlieb and Turkel 27 and 
Bayliss and Maestrello 28 . 
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The 2-4 scheme has been used successfully on a wide range of fluid and aeroacoustics 
problems 29 32 . Sankar, Reddy, and Hariharan 33 have evaluated this scheme for 
aeroacoustics applications. The solution procedure for the new scheme is as follows: 

In the present code, the optimized 4-6 Runge-Kutta time marching scheme of Hu, et. al 34 , 
is employed: 

Q k +4 _ l FBFB iBFBFBF l BFBF i_FBFBFBQk ^ 13) 


where k denotes the time level, L 4 and L6 denote four- and six-stage Runge-Kutta time 
integrations, and the superscripts B and F denote backward and forward spatial 
differencing. 


Note that the number of stages in the time marching procedure alternates between four and 
six. This procedure can be written as: 
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where the values of the upper coefficients are used in the four stage step and those of the 
lower coefficients are used in the six-stage step. Each derivative uses biased differencing, 
either forward or backward, providing inherent dissipation for the solver. Unlike the 
earlier MacCormack-type schemes, the stencil is not fully one-sided. 

Using a radial derivative at point j as an example, 


Forward: 


dQ 

dr 


-1 
A r 


0. 30874 Qj _! + 0. 63260/ - 1. 23300, +1 
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Backward: 


dQ 
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/ 0.30874Q /V | +0.6326Q, -1.2330Q/ V 

dr 

. _ A r 

i 
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(16) 


The sweep directions are reversed between each stage of the time marching scheme to avoid 
biasing, and the first sweep direction in each time step is alternated as well.. At the 
computational boundaries, flux quantities outside the boundaries are needed to compute the 
spatial derivatives, and these are extrapolated using third-order extrapolation based on data 
from the interior of the domain. 
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Due to the extremely sharp gradients of the shear layer near the inflow boundary, additional 
dissipation was required to damp out high-frequency errors in the jet shear layer. To 
accomplish this, sixth-order dissipation was added in the shear layer as: 

f Qj - 3 — 6Q/_2 + 1 5 Q/_i — 20 Qj 


f(q) = f(q) + 0.03 


+ 1 5 Q/+-| — 6Q/+2 + Oj + 3 + Qj - 3 
SQ j-2 + 1 5Qy_-j — 20 Qy + 1 5C?y + -| 


(17) 


^-60^+2 +Oy+3 


Computational Grid 

The grid used for these cases was 341 (radial) x 461 (axial) points. In the radial direction, 
the minimum spacing was Ar/D = 0.01 at the centerline, smoothly stretching to 0.13 at r/D 
= 16. In the axial direction, the minimum spacing was Az/D = 0.04 at the z = 0 boundary, 
smoothly stretching to 0.13 at the z/D = 35 boundary. The maximum spacing corresponds 
to 10 points per wavelength, which is well within the accuracy range of this code. All three 
cases used the same computational grid. A Courant number of 1.25 was used for these 
computations. 

Boundary Conditions 

Special attention is given to the boundary treatment in order to avoid non-physical 
oscillations which can render the computed unsteady solution unacceptable. Several 
boundary treatments were considered. 35 - 36 The boundary treatments discussed below 
were found to be stable, essentially non-reflecting, and suitable for the present jet noise 
computations. 
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At the inflow boundary (z/D = 0), the radial boundary is split into hydrodynamic 
disturbance and radiation regimes, which are treated differently as outlined below. 



At the hydrodynamic inflow boundary (r/D < 2), a small random disturbance is introduced. 
In space, the disturbance has the form of a 3rd order polynomial in order to obtain smooth 
derivatives in all directions, with a semi-random time component: 

e i(co+0 u )t 
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(19) 


A random number generator is used to determine the phase of each component of each 
mode and frequency. In this preliminary work, a Strouhal number of 0.2 was used. This 
Strouhal number was set using the conditions at the exit plane of the reference jet for 
compatibility. 


To introduce the input disturbance into the flow field, the time derivatives of the 
disturbance are added to the computed flow variables at each time step: 



boundary 



computed 
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In the hydrodynamic disturbance regime (r/D < 2), the Thompson inflow boundary 
condition is used 37 - 38 . In the Thompson analysis, the axial operator is decomposed into 
five 1-D characteristics. At a subsonic inflow boundary, four of these characteristics are 
incoming, and are set to zero for a non-reflecting boundary condition, while the fifth 
characteristic is outgoing and is computed from the flow solution: 

R, = p t '+pc{v z ') t = 0 

R 2 = p t '+C 2 p t ’= 0 

R 3 =pc{v r ') t =0 (21) 

R 4 =pc{v e ') ( =0 

R 5 = p,'-pc(v 2 ') t = {p,'-po(v z '),) compuied 

The five characteristic equations are then solved together to obtain the time derivatives of 
the variables at the inflow boundary. For a supersonic inflow, all characteristics are 
incoming, and all are set to zero. In order to prevent spurious oscillations, the amplitude of 
the outgoing characteristic was smoothly reduced to zero near the sonic line. 

Due to the specified disturbance at the hydrodynamic inflow boundary, the Thompson 
inflow boundary condition exhibited a problem in which some disturbances were convected 
in a radial direction and remained on the boundary for the rest of the computation. To 
alleviate this, the mean radial velocity was set to zero on the inflow boundary, and 
smoothly raised to the proper value by z/D = 0.7. 


tarnation 


In the radiation regime (r/D > 2), the conventional acoustic radiation boundary condition 
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applies: 


<fc=- vm 
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and M is the local Mach number. The spatial derivatives which appear in Eq. (22) are 
evaluated in an identical manner as the inner flow derivatives. 


Outflow Boundary Conditions 


The outflow boundary treatment is based on the asymptotic analysis of the linearized 
equations as given by Tam and Webb-^. The pressure condition is the same as that 
obtained by Bayliss and Turkel 40 , Enquest and Majda 41 , and Hariharan and Hagstrom 42 , 
namely: 


p't = -V(6) 


z , r , 
R P * + R Pr+ 



(25) 


However, for updating the rest of the primitive variables, Tam and Webb have shown that 
the momentum and continuity equations should be used to account for the presence of 
entropy and vorticity waves at the outflow boundary. The spatial differencing used in the 
inner code is employed to evaluate the derivatives which appear in Eq. (25). 
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For the outflow boundary at large radius and a local Mach number of less than 0.01, the 
outflow condition is replaced by the radiation condition described above. 

It must be noted that the Tam and Webb outflow boundary condition is formulated with an 
assumption that the mean flow is uniform, which is not true for the jet outflow boundary. 
However, the results given by this boundary condition are quite good, with very little 
reflection. 

Outer Radial Boundary Condition 

At the outer radial boundary (r = r m ax, 0 < z < z max ), the radiation boundary condition 
described above is used. 


Centerline Treatment 


In this code, the centerline boundary is represented with a point at the centerline, and a 
ghost point reflected across the centerline in the radial direction. Without azimuthal mode 
decomposition, the centerline treatment for a three-dimensional problem is not 
straightforward, and was addressed by Shih, et. al. 43 However, using the azimuthal mode 
decomposition method, the centerline boundary condition becomes straightforward: 
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First, Dahl's PNS solver was run to obtain the mean flows for the three jet test cases. 
These mean flows were interpolated to fit the LEE grid. The cases tested used a Strouhal 
number of 0.2 and modes n=0 and n=l for each jet. The input disturbance was identical 
for all three jets. In the interest of brevity, only the results of the n=l mode will be 
presented. 

Figures 3a-c show instantaneous snapshots of the pressure for the n=l mode of all three 
jets. Notice the preferred direction of radiation, and that there are two families of waves 
present, radiating at angles very close to the approximate angles given in Table 2. This 
shows the ability of the linearized Euler equation solver to capture both types of instability 
without any special treatment. 

Figures 4a-c show the contours of the radiated sound in dB for the n=l mode of all three 
cases. The magnitudes given are for comparison only, since the true magnitude of sound 
depends on the unknown amplitude and distribution of the inflow disturbance. Each 
contour represents a 5 dB change in the sound pressure level. It can be seen that the NVP 
jet disturbance reaches its peak farthest upstream, while the IVP jet disturbance continues to 
grow and peaks much farther downstream. 

Figure 5 shows the magnitude of the radiated sound as a function of the angle from the jet 
centerline. These results are taken at the outer boundary of the computational domain; thus, 
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the radius is varying from a minimum of R = 16 D to a maximum of R = 38.5 D. It is 
shown that the NVP jet is generally quieter than the reference single jet, while the IVP jet is 
as much as 8 dB louder. 


Conclusions 

This work is designed to illustrate the use of the linearized Euler equations for a parametric 
study of jet noise, as well as demonstrating the capabilities of the new MacCormack-type 
solver. The ability of linearized Euler equation solvers to use a mean flow from an existing 
CFD code in order to compute the noise radiated by instability waves in a supersonic jet is 
also demonstrated. 
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Figure 1 . — Mean flow velocity profiles, (a) For normal velocity profile coannular jet. (b) For 
inverted velocity profile coannular jet. 
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Figure 3. — Concluded, (b) In normal velocity profile jet. (c) In inverted 
velocity profile jet. 
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Figure 4. — Sound pressure level contours, (a) For reference jet. (b) For 
normal velocity profile jet. (c) For inverted velocity profile jet. 
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Figure 4. — Concluded, (c) For inverted velocity profile jet. 
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Figure 5. — Noise directivity results for coannular jet test cases. 
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